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Abstract 

We revisit the classification of polarization observables of the cos- 
mic microwave background. There exists a unified approach to the 
3x3 density matrix by which intensity, linear and circular polariza- 
tion are treated on an equal footing. The representations of the rota- 
tion group on the sphere contain certain right-left setup symmetries 
which have not been taken into account before. Left-right symmetries 
revise the construction of invariants and certain predictions based on 
symmetries. There are two true rotational invariants mode by mode 
in polarization data. Independent of models of cosmological pertur- 
bations, we emphasize methods to construct invariant distributions in 
order to test predictions of isotropy, axial or parity symmetries, and 
generate covariant and invariant statistics order by order. 



1 Introduction 

Current interest in the cosmic microwave background focuses on test- 
ing models of perturbations about an isotropic thermal default. Po- 
larization data is highlighted for providing new information about the 
relative size of scalar versus tensor classes of such perturbations. The 
groundwork for this task was set up years ago. And yet, surprisingly, 
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there are gaps in the current approach. For example, the distribution 
of CMB polarizations in various alternative cosmologies might be put 
into correspondence with the symmetries of the cosmology. Or one 
might try to assess deviations from uncorrelated isotropy using em- 
pirical distributions. Quadratic correlations do not encompass such 
questions. 

Here we revisit the kinematic task of describing polarization and 
the distribution of polarization, setting aside models of perturbations. 
The statistics of linear polarization are more subtle than commonly 
believed. The angular orientation of the polarization ip is measured 
at polar coordinates 6, 4>. Other variables that can be measured are 
the intensity (commonly called "temperature"), circular polarization, 
and the degree of polarization. Although the different variables have 
different transformation properties, we set up a framework in which 
all can be treated together under one unified description. 

The observables of polarized light are the wavenumber k = kk 
and polarization density (coherency) matrix pij =< EiEj* > for 3- 
component electric field E^. Since the field is transverse, it is conven- 
tional to reduce it to the transverse components, pjj =< E^E* 1 - > . 
This step complicates the full classification, and we do not use it. Next 
the astronomical conventions 1 employ unit vectors 6, (f> in the tangent 
plane to the sky, introducing a natural basis but also a troublesome 
coordinate-dependence. The 2x2 Hermitian matrix is expanded in a 
local complete set of tensors, via Stokes parameters: 



Here S^j is the 2x2 unit matrix orthogonal to the direction of obser- 
vation r. By construction here, 



tfj — OiQj ~ 4>i4>j'i 

t% = kh + 9jh 

tjj = 9i4>j - 4>i6j. 



(1) 



Q 



1 tr(pt E ) 



U 



1 tr(pt B ) 



(2) 



1 Actually astronomers use RA and DEC. 
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where tr indicates the trace. Born and Wolf use sq = I, s = Sj = 
Q,U,V. A pure polarization state oriented at angle ip relative to <p 
has Q = cos2ip, U = sin2ip. Under a rotation by angle 5 about the 
line of observation f, Q and U mix: 



The transformation property via twice the angle represents the tensor 
nature of the density matrix. 

Given some convention for interesting modes, one might predict 
them in a model for perturbations in the 9, 4> frame. The observer can 
then measure the modes in just the same coordinate frame to confront 
the model. Both model and observations then contain coordinate de- 
pendence, which stand to complicate or invalidate model-independent 
tests and tests based on symmetry. Towards dealing with these prob- 
lems, Zaldarriaga and Seljak (ZS) 3 and Kosowsky, Kamionkowsy and 
Stebbins (KKS)^2] developed elegant rules for projecting p E and p B 
into scalar fields. Recognizing that the linear polarizations in p 1 - are 
linear combinations of helicity ±2, ZS used the technology of New- 
man and Penrose IE] called "spin- lowering" to cast the observables 
into invariants. KKS project the 2x2 density matrix into a fiducial 
basis more directly, and show their method is equivalent. Our descrip- 
tion avoids spin-lowering and the entire 2x2 approach for technical 
reasons we will review. 

1.0.1 An Instructive Paradox 

The need for an appropriate framework for statistics of angular vari- 
ables is illustrated by an instructive paradox from biology 1 . An 
observer studying migrations records angles of bird flights 8i relative 
to East. The distribution is actually isotropic, as (say) 360 birds fly in 
360 different directions separated by one degree. The observer makes a 
distribution f(9) = dN/d6 from the data. Classical statistics suggests 
the average angle < 9 >= 9i/N is informative. For the isotropic 
distribution spread over < 9 < 360° he finds < 9 >= 180°, from 
which the observer reports that the average bird flies West. 

An obvious error is that 9 is a convention-dependent variable with 
the unpleasant transformation property 9 — > 9+5 when the coordinate 
system is rotated by 5. This problem is addressed by making explicit 




(3) 



(4) 
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the coordinate origin with an explicit angle 6 East- Then 9 — 9 East is 
invariant under change of origin. All higher moments of (9 — 9East) A 
are invariant. Yet despite invariance, such moments generally con- 
tinue to depend upon the convention for 0East^ The reader can ver- 
ify this with a few numerical integrals in the prototype distribution 
exp(—kcos(9 — 9 East))- F° r k = 1 the values (9 — 9East) 2 go from 
about 5.2 to 11.5, depending on 6$. Here 9q is an arbitrary parame- 
ter defining the zero of the coordinate. We call moments of this kind 
"convention dependent invariants" (cdi): they typify the troubles of 
circular statistics. 

Towards correcting such diseases one first decides upon a sym- 
metry group of transformations that will be used. One then maps 
the data into representations. Group representations are complete 
set of elements that transform linearly. Coordinate changes and the 
construction of meaningful invariants are greatly expedited. Yet the 
bird-flight example reminds us that invariance is necessary but not 
always sufficient for data analysis to be independent of conventions. 
Non-Abelian groups, in particular, have coordinate conventions that 
depend on ordering. 

For polarizations the 3x3 density matrix at each fixed position r 
is a reducible tensor of rank 2 under the spin rotation group, which 
decomposes 1x1^2+1 + 0, where the integers label the spin. Spin 
representations are labeled by a second quantum number s, thus 

S s 

ps, s = pij, 

where is a Clebsch matrix that depends on frame conventions. 
Then ps, s {8, <t>) over the dome of the sky has an expansion 

ps, s (9, 0) = 4^ m Yi(9, 0), (5) 

where Ym are the spherical harmonics. By construction these Ps sm 
elements transform like a tensor under a SO sp i n (3) x SO or bit a i (3) group 
of 6 parameters acting on the indices m and s, respectively. 

The actual continuous transformation group for rotations is SO to tai{3), 
whereby the 6 possible parameters of SO sp i n (3) x SOorutaii'S) are re- 
duced to 3 common parameters used in one joint rotation. There 
remains 3 fixed parameters associated with the ordering of opera- 
tions. There is hidden freedom in how the groups are joined together. 
Extra work, described below, is needed to cast Pij(9, 4>) into repre- 
sentations of SOtotai(3) taking into account all parameters. Our task 



4 



is to explain and explore the consequences of the extra freedoms for 
observables versus cdi. 

There is an interesting parallel with Wigner rotations known in 
defining spin-states within representations of the Lorentz group. More 
than one Wigner rotation can be used to make a basis! After a good 
basis is built, there comes a step of classification of its elements un- 
der other transformations such as parity. Good classification reveals 
convention dependence among invariants, which will otherwise hide if 
standards end at mere invariance. One can then categorize observ- 
ables hinging on invariant comparisons of two or more measurements. 
Finally, classification of the observables leads to classification of the 
distributions of the observables which allows symmetries to be tested. 

Our focus, then, is on making invariant distributions and using 
them to probe the symmetries of the CMB in a model-independent 
way. We pay careful attention to exposing human conventions that 
cannot be a feature of Nature. Yet we remind the reader in advance of 
basic validity in setting up coordinate systems full of conventions and 
comparing models with data in just the same conventions. Many pos- 
sible frameworks may serve so long as conventions drop out between 
theory and data analysis. The reverse is not true: not all "scalars" 
constructed in the existing framework are convention- free observables. 
For this reason the distinction between an "invariant" and an "observ- 
able" makes a difference depending on how data is analyzed. 

2 Transformation Properties 

We start with the 3x3 density matrix p(x) because it is what is 
measured. Let us review the transformation properties of such fields 
under SO to tai(3)- 

Let &a(x) be a set of fields transforming under an irreducible rep- 
resentation of spin S of the rotation group. With a rotation 
matrix, the transformation rule is 



^ s a (x)^^(x') = R s a ^$(x). 



(6) 



Let 




(7) 
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where 9 k are the parameters, a^g is the rotation generator for the spin 
of the field <J> a , and €ijk is the rotation generator on 3- vectors. The a 
matrices are matrix elements of rotation operators S, a^p = (a\ S k 
By definition 

[a\ a j } = ie ijk a k , (8) 

for i, j, k, cyclic. 

The infinitesimal change 5& a (x r ) = <& a {x') — <& a {x) is 

5<$> a (x) = 5<S> s a pin + 5<S>° a rbital ; 

^orUtal = -O^x'i-^jQvtf). (9) 



"dx>. 



This yields the orbital generator 



k _ . i d 

J-> — —l€ijkXi-^j, 



with the algebra 



[L\ V ] = ie ijk L k , (10) 



for i,j, k, cyclic. The SO sp i n (3)xSO or bi ta i(3) have now been exhibited: 
the final step is to unite them under SO to tai (3) defined by the diagonal 
subgroup transforming with the same parameters on both groups: 

J = L + S (11) 

The point of this review is that rotations employ two a-priori indepen- 
dent transformations accounting for the change in the spatial location 
(orbital) and the re-orientation of the indices (spin) implied by a ro- 
tation of a tensor field. 



2.1 Frame Fixing 

Since the SO sp i n {2>) and SO or bitai(3) are independent, there is nothing 
in mathematics to orient the basis of one relative to another: they as 
independent as angular momentum and isospin up to this point. 

Phase Conventions and Ordering: In beginning quantum mechan- 
ics one diagonalizes S z — > s z and L z — ► m. Choosing the same spin 
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Figure 1: A simple illustration of frame fixing. The coordinate lines for the 
location on the sphere (orbital angular momentum) have no a-priori relation to 
the abstract frame for spin-coordinates of "sticks". Although the convention on 
the right is not standard, it is as valid as any other. 

and space frame is a free convention since [L,S ] = 0. There remains 
a need to fix the phase of the eigenstates, which can be interpreted 
as the orientation of the x, y axes on the spin-world relative to the 
x', y' axes on the orbital world. This small point can be found in 
the book by Wiener [TT]: it creates a phase convention on each eigen- 
state. Yet since the eigenstates of J z = L z + S z already have phase 
conventions, and the rules of quantum mechanics make the conven- 
tions unobservable, the relative orientation of space and spin frames 
is invariably ignored. But the analogous problem on the sphere will 
become a central issue. 

For CMB polarizations it is more natural to diagonalize the helicity 
of the density matrix, the operator f ■ S. Notice however that f ■ S 
does not commute with L z : 

[L\S-f]=ie l3k ffr k . (12) 

Thus if we operate with f ■ S, followed by a rotation of the space frame, 
it is not the same as rotating the space frame followed by operating 
with r ■ S. Nevertheless S ■ r is a global scalar under the total rotation 
group, 

[J\S-f} = 0. (13) 

This causes eigenstates of f ■ S to be convention-dependent. In partic- 
ular, their phases depend on the convention and the order by which 
they are constructed. 
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To appreciate this dilemma, consider two natural procedures to 
set up coordinates on the sphere. Procedure 1 starts at the top of the 
sphere, the z axis, and diagonalizes L z and r ■ S = S z to make states 
\m, s) 1 . Let the usual x and y be orthogonal unit vectors to complete 
the definition of the spin basis. Any relative rotation about to the z 
axis of the space-and spin-frames leads to a phase which is ignorable. 
Define other states \ip,9,(fr} 1 by a rotation, 



9 j = exp(—iJ ■ 9) \m, s) 1 



Under the same rotation 



(14) 



S ^S' = RS 
r — > f' = Rf, 
S-f^S-f. 



(15) 



Thus 



S ■ f 



(16) 



With exp(—iJ ■ 9) = R(ip)R(9)R(cj)) this procedure is equivalent to 
parallel transporting the x, y, z basis to the conventional 9, 4>, f basis 
with a final rotation about the f axis by ip. 

Compare Procedure 2, under which we first rotate the spin basis 
by arbitrary Euler angles a, (3, 7, thus S l — > R(a, (3, 7)^ . Since S ■ f 
is a scalar, it is equivalent to rotating r % — ► r(a,@, 7) = R(a, 7)^^ . 
The frame defining orbital states is unchanged. Diagonalize f' ■ S to 
make states |m, s} 2 . Define other states \ip,9,(p) 2 by a rotation, 



exp(—i.J • 9) \m, s} 2 



(17) 



Under the same rotation 



S -» S' = RS 
f' _» f" = i?f', 
S-f'^S'-f". 



(18) 



Thus 



d • r 



(19) 
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This shows there are an infinite number of distinct ways to set up 
"a basis" while the helicity quantum number remains invariant under 
spin-rotations. One cannot ask Nature to know which procedure will 
be used. 

A little thought shows that there is no absolute criterion from 
mathematics alone to fix such conventions on a sphere; Figuren s h° ws 
a cartoon illustrating the problem. In much the same way, helicity 
states in relativistic field theory must be defined by some particular 
sequence of operations known as Wigner rotations. Starting from a 
particle at rest, a convention must be chosen to define its spin frame. 
This is followed by a rotation to orient the helicity states, which are 
finally followed by a boost. Since boosts and rotations do not com- 
mute there is an obvious convention dependence in the basis. Just 
as with SOtotai(3) being made from the direct product of two com- 
muting groups, the Lorentz group can be decomposed into commuting 
SU(2) x SU(2) to make representations. The relative orientation of 
the separate product group frames is an arbitrary convention, and 
nothing physical can depend on it. 



2.2 Left Cg) Right Setup Symmetry 

To build up representations it suffices to take outer products \s,S) (m, L\ 
in all possible combinations. We used freedom to chose one frame in 
the dual space, because it shows immediately that there exists unitary 
representations of the form U sm , UU< = U'U = 1, with S = L = j. 
These are the representations made famous by the 50(3) left-right 
sigma model. They transform under left-and right-groups as 

U - VlUVI 

where V L V ] L = V R V ] R = 1. Note that V L ^ V R in general. The 
freedom of spin-versus-orbital frame fixing amount to a Left-® Right 
setup symmetry in defining a basis. 

It is well known that unitary (j, j) representations of SO(3) make 
a complete set for function of three Euler angles. Defining as usual 

Di m (a, Pi) = {smj | U(a, /3 7 )) , (20) 

the completeness and orthogonality relations are: 

£ Di m (a, M)Di* m (a', 1 V) = ^jS(a - a')5(cos(3 - cos0)8{n - V), 

smj 
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/ 



dtlDl m {a, Pj)E%* m ,{a, A7) = g * ^ mm'Sss'Sjf- 



Here dQ = dt/jd cos 9d(f). Due to Left-Right setup symmetry it is an 
equally good expansion if we re-orient the spin-frames relative to the 
space- frames by an arbitrary rotation R 3 ss , (a A) > 7o) : 

Di m &, e, <p) = RiMPl ii)DLty, d > <t>)\ 

(22) 

We could have also rotated the orbital m indices, of course. There is 
one free-rotation for each j, so the freedom is considerable. 



2.2.1 Criteria for Observables 

For reference we define a true scalar 4>q as a quantity invariant under a 
one-time Left-Right setup SO (3) R(ao Ab 7o)j followed by continuous 
joint SO (3) R{af3, 7): 

(p Q -> R(a/3, i)R(a Po, 7o)^o-R t (" 0, 7) = <t>o- (23) 

This remnant of the original SO (3) x SO (3) transformations reduces 
6 continuous parameters to three continuous plus three fixed param- 
eters. 

We define a dependent scalar as a quantity unchanged when cxq = 
A) = 7o = are fixed. Observables that do not depend on the setup 
of spin and space frames must be true scalars. Dependent scalars can 
serve for convention-dependent invariants. 

The distinction of true and dependent scalars is a general feature of 
making product group representations. We can illustrate it intuitively 
with the group 50(2) x 50(2) (Fig. EJ). Let 4> be the azimuthal 
angle of a point on a circle in a 2-dimensional plane. Let % De the 
angle orienting a unit-length "stick" glued to each point on the circle 
- equivalent to a linear polarization. A continuous function /(</>, x) 
has the conventional expansion 

M X)=E fmne^e 2 ^. (24) 

mn 

Each exponential factor is an eigenstate of rotations about the origin 
in each product space. Since there are two independent origins of the 
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Figure 2: Frame-fixing for a direct product of two Abelian groups. The co- 
ordinate frames for position angles 4> versus "stick angle" 9 are independent. 
The facts of Abelian groups (see text) allows the issue to be ignored and the 
naive power to be a true scalar. 
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angles, the more explicit expansion including information on angular 
origins is 

f(4>, 00! X, Xo) = £ f mn e im ^-Me 2in ^\ (25) 

mn 

A change in the offset of either origin is the transformation 

f _> f p -™(0O-0o)-2m(xo-Xo) 

Jmn ' Jmn^ ■ v iL V 

Happily the usual power |/r?in| is ^ true scalar. In this case it does not 
matter whether one attends to the difference of partial or true scalars. 
Unfortunately this convenience does not extend to non-Abelian groups, 
where the order of operations determines the properties of the basis 
elements. 

Yet there is no need to overcomplicate matters. If one sought 
general covariance under arbitrary transformations, a coordinate grid 
where polarizations are measured might be an arbitrary curvilinear 
system independent of a curvilinear system locating the object posi- 
tions. But we are not discussing any transformations more compli- 
cated than rotations. It is enough (but necessary!) to deal with the 
complication of frame fixing. In the next Section we discuss how the 
convention dependence originated. 

2.2.2 Parity Versus CMB s-Flip 

Classification is completed by parity P. The conventional parity trans- 
formation for a tensor operator t(x) of spin S is 

t(x) -> Pt(x)P t = (-lft(-x). 

Thus P = PsP(L) where P = Ps and P = Pl are the spin- and space- 
parity, respectively. Angular momentum J commutes with parity: 

PJP ] = J; (27) 
[P,J]=0. 

Hence parity is a scalar operator. Position anti-commutes: 

Pfpt = (28) 
{P, x} = 0. 

It is elementary but not uncommon to see the signs of Eg 1271 1281 re- 
versed, which makes a difference in the following. 
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True parity transformations of CMB density matrices relate the 
data on opposite sides of the sphere. In the CMB literature the oper- 
ation of reversing S ■ r with fixed direction f is often called "parity," 
but since r reverses under parity (Eq. 05|) we prefer "CMB s-flip" ir s . 
Then 

ir s S ■ ri: s = —S ■ r; (29) 
ir s L ■ rir s = +L ■ f, (30) 

The algebra implies that quantum numbers s — > — s, m —* m: 

S ■ fir s \s, m) = —tt s S ■ r \s, m) = —stt s \s, m) . (31) 

One defines E and B modes as even or odd under s-flip: 

D% = -{D{ m + Di 2m )/2- 

tt s D e / b = ±D E ' B , (32) 

where the slash indicates the E or B case. 

Now as a consequence of Eqs. I29| 1301 the CMB s-flip is not a 
scalar under rotations. Since the E and B modes are made from spin- 
frames of a particular convention, they retain the convention. This 
is intuitively apparenty if one temporarily forgets about the global 
pattern of E and B mode "cartoons," and simply looks along a line 
of sight in some fixed f direction. Locally one sees a single stick 
whose orientation is pure convention relative to an arbitrary local 
coordinate system. The data is local, but the classification is not! 
The transformation of E and B modes into one another under spin- 
rotations is simply Eq. Elonce again. 

There are several reasons for missing this simple transformation. 
First, true parity is a true scalar operator that is very familiar, and 
classification by parity cannot change under full rotations. By re- 
placing parity with ir s one tends to take true invariance for granted. 
Second, conversion of D sm — > D^, eliminates the s index which 
transforms, creating a notation with appearance of true invariance. It 
is illusory because the phase of the S ■ f states is free, which is just 
the same phase that varies under under spin-rotations, and just the 
same phase fixed by hidden conventions. Earlier we mentioned that 
— S-r is not a true scalar under either the spin or orbital groups, 
but only jointly. Classifying E and B using spin-parity relative to 
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the spin-group causes tt s not to be invariant under Left-Right setup 
symmetry. 

Finally it is easy to confuse CMB s-flip with spin-parity P s and or- 
bital parity P^. There exists a notion of "CMB m-flip" 7r m completing 
the algebra: 



Note the signs. As with eq. |^ this algebra implies that quantum 
numbers s — > s, m — ► —m. Overall parity P = PsPl = K m ir s , by 
which the product makes the usual scalar operator. 

Why are E modes and B modes nonetheless thought to be invariant 
classes? It is because they are dependent invariant classes, invariant 
only under fixed spin-frame and orbital-frame transformations, not 
fully invariant classes. 

Why don't these intricacies plague conventional classifications? It 
is because spin-parity and orbital parity are true scalars, not referring 
to a spin-frame convention. 



Here we review background related to the frames and their history. 

3.1 Euler Angles and D Matrices 

The simple and natural coordinate system used in astronomy 2 happens 
to correspond beautifully with the formal procedures used in defining 
our basis. 

Take the XYZ coordinate frame with origin O and translate it 
without rotation to the top of a unit sphere. Rotate by angle <fi about 
OZ, so that the new unit vector x is locally parallel to 9. Rotate by 
angle 9 around Y' in this system, sending the axis Z — ► Z' located 
along 9, <j). Rotate around the new Z' axis by ip until it is parallel 
to the linear polarization. The sequence of rotations is denoted by 
the operator R((f), 9, ip). It is easily shownjl] to be equivalent to first 
rotating by ijj about OZ, then 9 about OY, then <f> about OZ. This 

2 Sign conventions of course must be taken into account 



n m S ■ fir m = +S • f; 
tt s L ■ fir s = —L ■ f. 



(33) 
(34) 



3 Background Review 
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allows identification of the matrix elements, and in particular diago- 
nalization about the fixed axis OZ, with the result 

DU^ 6, VO = e^4(fl)e^. 

The di m matrix elements are attributed to Wigner, who developed an 
explicit formula: 

Aj + m)\(j - m)\ 1/2 
1 {j + s )\{j-s)\ J 

Sp, l)3 +s-r ( 3+S \( j-S \ 

Z_^ v y r Jyr — s — mj 

e is ^sin 2j (9/2)cot 2r ~ m ~ s (9/2)e im ^. 

For practical data analysis it is straightforward to enter the D° sm 
formula into a computer once, transfer the angles 4>, 8, ip for each data 
point, and generate a corresponding array of representations for each 
j- 



Di m {4>,6^) = 
x 
x 



3.1.1 Spin-Lowering 

In 1959 Jacob and WickfT^j introduced the helicity formalism for the 
representations of the Poincare group. Early papers on helicity con- 
tained detailed discussions of the Wigner rotations and phase con- 
ventions. In 1966 Penrose and NewmanjS] (NP) applied the helicity 
formalism to a problem in gravitational radiation. For unknown rea- 
sons reference to Jacob and Wick, and any discussion of frame-fixing 
were omitted. Instead a notion of "spin-lowering" was introduced. 

Spin-lowering is the addition of units of orbital angular momentum 
to the spin angular momentum to make objects of net helicity zero. 
The NP operators 8 operating on an object of helicity s are defined 
by 

3 = sin s e(— + icscO— )sin' s 9. (35) 
oO ocp 

Lowering is done by 3L The interior bracket is simply the ordinary 
orbital angular momentum raising and lowering operators L x ± iL y 
relative to a basis of x = 6„ (j) = y. Thus spin lowering (raising) 
subtracts (adds) a unit of f(9, 4>)-L. The factors oisin s are Jacobians 
needed for operating on objects of helicity s (called "spin-weight" by 
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NP). As above the orientation of the (x,y) or 6, <p spin-frame is 
a convention. NP remarked cryptically that "9 is invariant under 
transformations that preserve the metric on the sphere". This is the 
particular phrase that hides the problem: The metric of humanly- 
defined frames has no pretense to invar iance, and is not an observable. 

./VP's program continued with objects constructed by lowering or 
raising the spin-weight to zero and expansion in spherical harmonics. 
Integration by parts clarifies this: 

{lm | 9 s p 8 ) = /s st Im | p s \ . (36) 

The left hand side shows the raised or lowered states expanded in 
states of orbital angular momentum I, m. The right hand side is the 
expansion of the original object in a basis defined by operating with 
9 s t. The basis of |S s ^m) were called the spin-weighted harmonics 

s 1 ra- 
it seems the s Yj^ initially appeared to be something new and sig- 
nificant. However Goldberg et al6\ soon showed that the expansion 
was simply an ordinary tensor spherical harmonics expansion in the 
helicity basis. The identity needed is simply 

27-1-1 1 / 2 

s Yl(8, <f>) = Dl sm {4>, 9,ip = 0). (37) 

Notice that a convention-dependent zero-point of ipo = appears. 
This is an SO(2) tip of the iceberg of frame-fixing. By omitting any 
mention of the angle t/j and the need for frame-fixing, NP seem to 
have created a misleading impression, to whit: one can misinterpret 
Eq. as projection onto an absolutely meaningful basis free of con- 
ventions. Yet there exists an infinite number of equivalent bases. The 
distinction of spin and orbital frames also underlies our use of index j 
for total angular momentum, as opposed to the inappropriate symbol 
I for orbital angular momentum in some absolute frame. 

Electric and Magnetic Modes: A complementary construction of 
tensor spherical harmonics was given by Zerillii in 1970. Zerilli[7] re- 
introduced the terminology of "electric" and "magnetic" linear com- 
binations that had been coined by Matthews [8]. Yet Zerilli failed to 
emphasize that the parity of the basis elements is simply a mathemat- 
ical detail having nothing to do with the parity of objects expanded 
within it. ThornejS], for example, used both the electric and magnetic 
basis elements to expand gravitational radiation for a metric of even 



16 



parity. Similarly, in electrodynamics |1U| the vector spherical harmon- 
ics are confusingly named electric and magnetic while a general electric 
field (say) needs to be expanded in both electric and magnetic modes. 
It should be clear that the conventions of "electric" and "magnetic" 
basis elements represent nothing observable by themselves, and the 
freedom of frame-fixing of spin and orbital bases remains unchanged 
no matter how the elements are re-arranged. 

3.1.2 Gradients and Curls 

KKS avoid spin-lowering by contracting pij with conventions for basis 
elements made from the derivatives of Y^. They demonstrate that 
their results are just the same as ZS due to an integration by parts 
(Eq. I36|) . The hidden frame- fixing convention lies in using the deriva- 
tives defined via the oribital frame to make the spin frame, overlooking 
the Left-Right setup freedom. Obviously many other conventions ex- 
ist, and the evident correlation of the spin- frames and orbital frames 
(oriented along the conventional polar axis) chooses a convention for 
Left-Right setup symmetry. 

The method of KKS makes the arbitrariness of frame fixing ex- 
plicit, while it is complete hidden in the NP method, generating a 
small paradox. To support belief that the basis has an absolute mean- 
ing, KKS suggested the combinations of derivatives V and f x V should 
be viewed as decomposing the field into "gradient" (G) and "curl" (C) 
modes. The authors cite the Helmholz theorem that such a decompo- 
sition should be unique. 

However there is no parallel of the flat space Helmholz theorem 
on the sphere. The analogy fails even for the s = 1 vector spherical 
harmonics. Suppose v(x) is a vector field on open flat 3-space, and 
one seeks </>, A such that 



v = V</> + V x A. 
In flat space this is solved 



(38) 



V x v = -V 2 A + VV • A; 
V • v 

0=—^; A = v-V<j>. 



(39) 
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Such Green functions have infinitely many symmetries 1/V 2 — > 1/V 2 + 
h, where V 2 /i = determines homogeneous solutions fixed by bound- 
ary conditions. The usual Coulomb kernel is then fixed by choosing 
solutions that vanish at infinity. Yet there is no analog of "infinity" 
on the sphere, and all solutions are doubly periodic. Hence there is no 
unique definition of 1/V 2 -, and every effort to construct one creates a 
convention, which is the frame-fixing in yet another form. 

We construct and discuss the Green function 1/V 2 - in an Ap- 
pendix. For an intuitive understanding of the breakdown of unique- 
ness, recall the physical interpretation of 1/V 2 - is the potential of 
a point charge, with electrostatic field — Vr/V 2 -. Gauss' Law then 
requires the point charge be cancelled by another an equal and oppo- 
site charge somewhere on the sphere. It is simple, but by no means 
mandatory, to place a canceling point charge on the opposite side of 
the sphere. This creates an axial anisotropy convention... indeed a 
new axis is introduced for each and every source! The corresponding 
spin-frames are the unit vectors of the electric field from the initial 
charge to the arbitrary location of the mirror charge. When the arbi- 
trary mirror charge is re-located the arbitrary spin-frames are revised: 
and therein lies the frame-fixing. 

The failure to have an ordinary Helmholz theorem is topological 
and unavoidable. It is interesting that this problem has been noticed 
in oceanography ^3]- The authors of Ref. ^3] state "The decomposi- 
tion of an eddy flux into a divergent flux component and a rotational 
flux component is not unique in a bounded or singly periodic do- 
main... assertions made under the assumption of uniqueness, implicit 
or explicit, may be meaningless." We don't say the process is mean- 
ingless. Conventions can be fixed to make it meaningful. However 
convention dependence must be brought out when comparing to ob- 
servables. 

4 Classifications 

Here we review the classifications of the density matrix in our scheme, 
towards developing convention- free invariants. 

4.0.3 Classification of p 

Under the spin subgroup p ~ 1 <g> 1 is reducible, and >2 + l + 0, 

where the integers stand for S. Let (S, s) stand for spin and helicity. 
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Then Q±iU ~ (2, ±2), V ~ (1,0), I ~ (0,0). Other components are 
absent because p is transverse. Since we are using p as a 3 x 3 matrix, 
the wave number k is the eigenvector of p, oriented along — r. 

Under spin-parity the electric field's vector E — ► —E, so that p — > p 
is even. The spin-1 associated with V is a pseudovector because the 
Clebsch for 1 x 1 -> 1, namely e^, needs a handedness convention. 
From products of p and k ~ (1,0), and using the addition rules of 
angular momentum, there is one pseudoscalar, Stokes V = tijkPijkk- 
From p we can make one new scalar beyond Stokes I, the degree 
of polarization \J Q 2 + U 2 . This finishes the local classification of p. 

4.0.4 The Data Map 

Given a set of measured angles in data <f>i 9{ ipi , the empirical distri- 
bution is a sum of delta functions concentrated on data values. The 
empirical distribution fdata is expanded in the basis of D's as follows, 

f((j>i 9i %ln) = "^2~ Am D im\ 

i 

P i sm = Y, DJ s*m(^0i^). (40) 

i 

Procedures for "cleaning" or interpolating data may appear here to 
incorporate features of instrumentation or to improve statistical reli- 
ability. Whatever method is used, we call Eq. blithe data map from 
the angles to the coefficients. 

Frames from data: Interestingly, data summarized by pi m contains 
within itself certain preferred spin and orbital frames. Statistical anal- 
ysis of the different frames may reveal new insights into the CMB, just 
as the orientation of the ordinary CMB multipoles has led to new in- 
formation. 

The "canonical frames" are found by diagonalizing pP sm for each I, 
writing 

pL = «j tt A^(^ T )« } (41) 
where T denotes the transpose. In coordinate-free notation 

pi = ^2\ u ja ) X ja (v ja \ . (42) 

a 
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This is the singular value decomposition (svd), which is unique. Here 
h a i = A^, A^ 2 ^ are the singular values ("eigenvalues"), invariant un- 
der SO sp i n (3) x 5 , O or ftj ta ;(3), and therefore invariant under the full ro- 
tation group. Matrices u J s a , Vm are the orthogonal frames constructed 
by diagonalizing the Hermitian combinations p> p*^ and pP^ (P: 

pVw = (^W^'W. (43) 

For CMB density matrices classified under s = 0,1,2 there are at 
most two orthogonal terms, Vm\ Vm for each j, while the range of 
m = 0, 1, 2, 3... is unrestricted. As eigenvalues of a positive Hermitian 
matrix, the diagonal elements (A J ) 2 are real and positive, and A? are 
conventionally defined to be positive to fix the phases of u{ a and Vm . 

The diagonal direct products show that any data for pi m can be 
interpreted for each j as two true invariants, plus inherent orienta- 
tion of spin frames strictly correlated with inherent oriented orbital 
frames. The orbital frames of course contain more information than 
the invariants: they are covariant under rotations. The invariants are 
readily extracted from ordinary traces of the matrices: 

(A (ib')2 + (A (2)^ = ^yt ); 
(A (i)j )2(A (2)j )2 = det (p? p jt) ) (44) 

where tr and det denote the trace and determinant. These are equiv- 
alent to the invariants exhibited earlier. 

4.0.5 Intrinsic Frames of Basis Elements: 

Just as above, the convention-dependence of the basis states can al- 
ways be discovered by the svd procedure. Given elements D{ m in a 
basis labeled 'V , there is a unique expansion 

DL {v) = < (*0A J 'V a ))» (")■ (45) 

The particular orthogonal matrices encode the net history of parallel 
transport used to construct the particular basis, which "latches onto" 
unit vectors (v ]a ) m and "returns" U$ clS 9i result. 

Hence if one experimenter generates a set of D 3 sm matrices accord- 
ing to a given convention, another experimenter can find the intrinsic 
frames, which are equivalent to the ordering and phase conventions. 
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5 Distributions 



We now turn to the construction of manifestly invariant distributions 
of the density matrix. An invariant distribution must be a function of 
an invariant argument / = f(p, n), where n represents tensor param- 
eters of appropriate type. When cast into our basis, / = /(psm> n J sm ), 
where v? sm are the matrix elements of parameter p. It is useful that 
each subspace of index j is independent. With distributions one can 
examine various correlations defined as moments, or use other tests 
such as likelihood, depending on the complexity of description sought. 
In order of decreasing symmetry, we list various cases useful to cate- 
gorize data. 

5.0.6 Isotropy and Parity 

Isotropy requires f(p, p) — > f(p) = f(I, Q 2 + U 2 , V), since these 
arguments are the only invariants made from p. If the distribution is 
invariant under parity, then f(I, Q 2 + V 2 , V) = f(I, Q 2 + U 2 - V), 
yielding < V >= 0, and similarly for odd powers of V. 

We have shown that E and B modes are not true invariant classes. 
There comes a paradox that parity symmetry has been used to make 
predictions of E and B mode correlations. Examine the bilinear cor- 
relation 

^Is', mm' =< PsmP's'm' > ' (^) 

where * indicates complex conjugation. Isotropy requires that C 3 SS , mm , 
be an invariant tensor. 

The only rotationally invariant tensor symmetric in ss' and mm', 
with generally different dimensions, is given by 

CL,, mm/ = CU ss r5 mm ,. (47) 

We then find the correlation 

< p% j p m ' j * >= CH mm , (5 22 - = 0. (48) 

That is, we find 

< ffttfftl* >= 0. (49) 

This a consequence of isotropy without mention of parity symmetry. 
In the current literature Eq. 0H1 is identified as a consequence of ir s 
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symmetry, not true parity. Perhaps 7r s -invariance might be something 
reasonable to propose, but it lacks the fundamental character of true 
parity. It does not break parity to find < Pw'Vro*'* >^ since the 
actual CMB is not isotropic. 

5.0.7 Axial Dependence 

The next simplest symmetry class would contain axial parameters 
n, which might be vectors (oriented) or the eigenvectors of tensors 
(unoriented). There is an orbital scalar made linearly via f-h. One can 
make a spin-pseudoscalar e-ijh^iPij which depends on Stokes V, and a 
quadratic spin-scalar by the combinations pijhihj. Continuing, from 
a single h there is an infinite number of tensors h<S)fi...h of arbitrary 
order, which can be reverted to objects transforming like helicity s or 
orbital angular momentum m: finally tensors n° srn in each subspace 
j. Current data for the CMB does have a symmetry-breaking axis, 
obtained from the dipole (j = 1) intensity map. Several authors [Tl] 
have noted that a nearly coincident axis appears for j = 2, 3, 4. Our 
analysis suggests that an axial dependence might be related to the 
< TE > correlation WMAP observes [IS]. 

When isotropy fails, Eq. 03 fails. It is clear that correlations of 
the type < p , amp 3 s * m i > are not manifestly invariant. The rotationally 
invariant quadratic correlations are of the form 

c j = Y,< pLpfa > • (50) 

sm 

The summands are linear combinations of the invariants listed earlier. 
By comparing this invariant with C 3 SS , m , various test of isotropy can 
be undertaken. 

Now suppose some signal violating isotropy is observed: how does 
one proceed? By analogy with familiar Gaussian correlations, the 
C 3 , , serve as model parameters one can use to build up heuristic 
distributions. Thus one can made simple likelihood tests, comparing 
/ = f(p) ( tne isotropic default) with / = f{p 3 smC 3 ss , mrn ,p l s * rn ,) 1 which 
uses an invariant built from the data itself. (We hasten to add that 
data sets are enormous and complex, and data analysis itself is an 
art.) 

We note that in the presence of anisotropy the relative orien- 
tation of the space and spin-frames can obviously be fixed in an 
absolutely meaningful way. In Cartesian language, we may define 
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Q n ~ hiPijhj] U n ~ fiiPijejpqhpkq. In this case Q and U can be cast 
into physically-based coordinates, not from p alone, but from p rel- 
ative to h. Now with U n a pseudoscalar, parity symmetry with the 
anisotropic distribution will predict < Q n U n > = 0, a meaningful re- 
lation for E and B modes defined relative to h. This explains why 
cartoons of E and B modes have visible parity properties: they are 
cartoons emphasizing a particular h axis. 

5.0.8 L ■ S Coupling 

The most familiar coupling of space- and spin-angular momentum 
in quantum mechanics is called the U L ■ S" coupling. The analogue 
for CMB data is straightforward to construct. Let S),J be matrix 
elements of the spin- angular momentum generators. Only the s, s' = 
±2 components will actually be used. Let L^J be matrix elements of 
the spin-j generators - these are not gradient operators, but matrices 
of numbers. Then the following is a dependent scalar for each j, 

c j LS = Y. tr ip i] s {j) -p il{j )> 

k 

= E &STPs'rn>L% k m (51) 
mm'ss'k 

This represents a correlation of the global polarization and angle de- 
pendence on the dome of the sky. The expected value is zero when 
there is isotropy in either the spin or space dependence. 



5.0.9 WMAP 

The observation of correlations that are convention-dependent is by 
no means useless information. Results so far reported by WMAP [15] 
include a "< TE > lr ' correlation showing surprising structure at small 
I. Sometimes the correlation is defined by 



C{9) 



where t-qi and r\- are location where Q and I are measured separated 
by angle 6, and symbols Wi, Wj are weights introduced to correct the 
data. The modified Stokes Q' is defined as the Stokes Q parameter 
relative to an origin oriented along great circles connecting pairs of 
points i, j 
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One can understand WMAP's procedure as choosing spin- frames 
at TQ t i via a fiducial z axis at r/j. Fix index j; the great circles are 
simply the coordinate lines of local Oj . With this convention calculate 
the unvarnished £7 modes of Eq. [2] relative to 9jj. Continue, moving 
to the next j and sum. Does this make a convention-free observable? 

Unfortunately no. From our analysis of frame-fixing there is one 
free choice of frame orientation to fix per I. This convention is a con- 
vention however it is done. For instance one might choose to measure 
Q' relative to the normals to the geodesic; then angle tj) 1 — > ir/2 + tp', 
whereupon < TE >'—► — < TE >\ We do not believe that the sign 
of a physical observable is an arbitrary thing I by I. More directly it 
is just arbitrary to develop the frames by parallel transport from fj j. 
(Suppose for instance that all the data came from a small region of 
the sky, as in earlier studies not covering the full sky. Then a bias 
relative to this choice of z is obvious.) We have shown no absolute 
frame fixing exists, and that < ET > correlations are not true scalars. 

As mentioned earlier, comparison of theory and data treated within 
any fixed and consistent procedure can always be made. WMAP in- 
terprets their correlation data as signals of a late ionizing phase, citing 
several theoretical predictions. We are unable to assess consistency of 
data analysis and theory which is liable to be quite complicated. Yet 
breakdown of true invariance gives no immediate reason to doubt the 
consistency of the comparison reported. It would certainly be inter- 
esting to explore the hidden freedom of convention when comparing 
data and models. 

It would also be very interesting to see more varied and versatile 
CMB polarization data analysis to complement more familiar power 
spectra. Testing distributions based on the invariant and covariant 
decompositions suggested earlier may have a rich potential for discov- 
ery. 
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6 Appendix 

In this Appendix we show analytically there exists no analog of the 
Helmholz theorem on the sphere. 

We may show this by contradiction. Consider a region "1" where 
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Vt-^ti 7^ 0, where T represents the transverse components in the tan- 
gent plane. Use Gauss' Law to integrate over the region, d 2 xx Vr ■ 
v T1 — > $ dSx-VTi- The same contour is the oppositely-oriented bound- 
ary of the complementary region "2" of the sphere. Hence 

j> dSr ■ vti- + <j> dSr • vt2 = 0, (52) 

meaning that the integrated divergence vanishes. 

We now seek an operator 1V^ such that we recover the sole diver- 
gence of region 1, 

This is a contradiction, because the net divergence must be zero, while 
the net divergence in region 1 is non-zero. The integration over region 
1 will instead assign a corresponding anti-divergence to region 2. What 
rule shall be used? Without actually integrating over region 2 (the 
region 1 integration probing only the boundary) it is impossible to 
recover the correct divergence coming from region 2. 

A specific example: consider the divergence singularity of v — > 9 
in the northern hemisphere. The only possible scheme to re-distribute 
this singularity in a rotationally symmetric way must involve a con- 
stant divergence. This defines a Green function Go such that 

V%G (6, 8', 0, 0') = 5(cos9 - cos9')5(<f> -(/>')-—■ 

The integral of the right-hand side is zero, as required, but we have 
an operator which differs everywhere from the desired Green function. 
It is clearly possible to force the canceling divergence to be an iso- 
lated singularity of opposite "charge," but where on the sphere shall 
it be placed? Thus all Green functions contain arbitrary symmetry- 
breaking conventions: introducing an artificial axial anisotropy is the 
most symmetric result possible. 

Construction of the Green functions goes as follows: Define G(f, f') 
such that on the sphere 

V 2 T G{f, f') = 5 2 {f, f'). (53) 

With f(9, 4>) in a reference coordinate system, expand G in spherical 
harmonics, whereby V\ — > 1(1 + 1) on each mode. Completeness 
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suggests the usual solution dividing out the operator, namely 

m,l=l v ; 

Note the I = term is divergent and had to be omitted. It follows 

V^G = 5(cos6 - cos6')5(<t> - <j>') - —, (55) 

which smears an anti-charge uniformly over the sphere. To compen- 
sate add any solution h such that V^/i = I /Air. Proceed by expansion, 

v 2 T h = ^i(i + x ) h im Yi(e, <l>)Y£(e', </>') = i-. (56) 

Im 

Note that either the expansion breaks rotational symmetry with some 
non-trivial tensor hi m , or h\ m — > hi does not depend on m. In the 
second case the addition theorem for spherical harmonics gives 

+ <t>)Y%{e',4>') = 

Im 

\- (21 + 1)1(1 + 1) hi . 1 

where 7 is the angle between f and r' and P; are the Legendre poly- 
nomials. Comparing the left and right side, 

7 ^0 



(2/ + + 



meaning that hi does not exist to satisfy the requirements. This leaves 
the first case, hi m 7^ 0. An infinite number of solutions exist but all 
break rotational symmetry. For example V^(Y^ — Yq) = for all / 
and any to. 

We can apply this calculation to the discussion of charge-and anti- 
charge of Section 1.21 The ansatz is 

^ Yl(t>, 4>)\±(t)', <Q Yj[* - e, -way, #) 

I = £„T^ W+i) • (58) 

whose Laplacian develops two 5 function singularities with the offen- 
sive 1/4-7T terms canceling. 
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